Quantum electrodynamics near a Huttner-Barnett dielectric 
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We build up a consistent theory of quantum electrodynamics in the presence of macroscopic po- 
larizable media. We use the Huttner-Barnett model of a dispersive and absorbing dielectric medium 
and formulate the theory in terms of interacting quantum fields. We integrate out the damped 
polaritons by using diagrammatic techniques and find an exact expression for the displacement 
field (photon) propagator in the presence of a dispersive and absorbing dielectric half-space. This 
opens a new route to traceable perturbative calculations of the same kind as in free-space quantum 
electrodynamics. As a worked-through example we consider the interaction of a neutral atom with 
a dispersive and absorbing dielectric half-space. For that we use the multipolar coupling /i, • D 
of the atomic dipole moment to the electromagnetic displacement field. We apply the newly de- 
veloped formalism to calculate the one-loop correction to the atomic electron propagator and find 
the energy-level shift and changes in the spontaneous decay rates for a neutral atom close to an 
absorptive dielectric mirror. 

PACS numbers: 31.70.-f, 41.20.Cv, 42.50.Pq 



I. INTRODUCTION. 

Quantum electrodynamics is a well-functioning the- 
ory which accurately predicts a wide range of phenom- 
ena, not just in high-energy physics but also in atomic 
physics. The best known quantum electrodynamic effect 
in atomic physics is certainly the Lamb shift which has 
by now been calculated to very high accuracy |l| . If the 
atom is located not in free space but instead near a re- 
flecting surface, which could be dielectric or conducting, 
then the reflection of photons from that surface leads 
to the Lamb shift acquiring a distance-dependent com- 
ponent, the Casimir-Polder shift, whose gradient yields 
the Casimir-Polder force between atom and surface. Al- 
ternatively, the Casimir-Polder shift can be viewed as a 
Stark effect where the role of the electric field is played 
by the non-zero and position-dependent electromagnetic 
vacuum fluctuations in the presence of dielectrics 0- 
In order to study the Casimir-Polder effect and related 
quantum electrodynamic effects due to the presence of 
macroscopic material boundaries, one needs a theory 
of the quantized electromagnetic field in the presence 
of such boundaries. The method of field quantization 
largely depends on how sophisticated a model of the ma- 
terial's optical response one assumes. In the simplest 
case one might assume perfect reflectivity of the surface. 
The quantization of the electromagnetic field can then 
be achieved quite easily by a normal-mode expansion of 
the field, where the electromagnetic field is expanded in 
terms of a complete set of solutions of the homogeneous 
Helmholtz equation. The presence of the boundaries is 
taken into account by imposing appropriate boundary 
conditions on the electromagnetic field. Quantization 
is then accomplished by promoting the expansion coef- 
ficients of each mode to creation and annihilation op- 
erators which are required to satisfy bosonic commuta- 
tion relations. This approach of canonical quantization 
has the advantage of being simple and therefore workable 



even for complex geometries Q but the perfect-reflector 
model for the surface lacks essential physical features, e.g. 
evanescent modes, which may have a dramatic effect on 
predicted quantities An improvement is to consider 
the material as a non-dispersive and non-absorbing di- 
electric characterized by a single real number, an index of 
refraction. Then field quantization can still be achieved 
by canonical quantization using field modes, although the 
specific implementation of the method requires a lot more 
care than for perfect reflectors Q. 

Canonical quantization of the electromagnetic field in 
terms of normal modes runs into difficulties when one 
wants to include in the formalism realistic properties of 
dielectrics. The response of the material's surface to the 
electromagnetic radiation in reality depends on the fre- 
quency of the impinging radiation. Furthermore, causal- 
ity requirements demand that any dispersion is always 
accompanied by absorption. However, a naive incorpora- 
tion of absorption into canonical field quantization leads 
to field commutators decaying in time, i.e. an inconsistent 
theory. Therefore, in any model of interaction between 
real dielectrics and the electromagnetic field, the field has 
to be coupled to a reservoir in order to simulate the ab- 
sorptive degrees of freedom 0] . This can be done in a 
number of ways. One is to model the absorptive degrees 
of freedom by adding to the operator-valued Maxwell 
equations Langevin-type fluctuating noise-currents that 
ensure that the canonical commutation relations do not 
decay in time but rather take the expected form Q. In 
this approach the field equations are solved by using the 
Green's function of the wave equation, and the noise- 
current operators and their properties play a major role 
in describing the dynamics of the coupled field-dielectric 
system. A number of papers have provided an a pos- 
teriori microscopic justification of such a procedure by 
deriving the commutative properties of the noise-current 
operators that were otherwise introduced ad hoc 

A more direct approach to modelling the interaction 
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between the electromagnetic field and an absorptive di- 
electric is to explicitly include from the outset in the 
Lagrangian (or Hamiltonian) the matter degrees of free- 
dom that are responsible for absorption. The dielectric 
is then envisaged as consisting of a continuum of har- 
monic oscillators coupled to a reservoir which consists of 
yet another set of harmonic oscillators. This quantum 
model of a classical dielectric was originally introduced 
by Hopfield [11] . The first Fano-type diagonalization 
of the resulting Hamiltonian was achieved for fields in 
three dimensions in Ref. for a bulk dielectric and the 
general treatment of inhomogeneous dielectrics followed 
in Ref. This model has also been extended to in- 

clude spatial dispersion [l^ and magnetodielectrics [iBl ■ 
Practical applications of the Huttner-Barnett model, e.g. 
the calculation of spontaneous decay rates , work well 
for bulk dielectrics where simple forms of the relevant op- 
erators can be found, though an additional difficulty is 
that in a bulk medium local field corrections play an im- 
portant role and need to be included. On the other hand, 
complications that arise due to inhomogeneities of the di- 
electric have previously led to unwieldy and impractical 
results; the conceptually very interesting work by Yeung 
and Gustafson [l^ uses the Wiener-Hopf method to cal- 
culate the photon propagator of the vector field A in the 
presence of an absorbing dielectric half-space, but the 
result is so complicated that it has to be Fourier trans- 
formed and evaluated numerically, whence all subsequent 
calculations are also necessarily only numerical. 



In this paper we demonstrate that by starting from a 
Power- Zicnau-Woolcy type of Hamiltonian rather than 
adopting minimal coupling, one can carry out explicit 
and easy-to-follow perturbative calculations in quantum 
electrodynamics in the presence of an inhomogeneous 
Huttner-Barnett dielectrics. We apply the formalism we 
develop to the problem of calculating the energy-level 
shifts and change in spontaneous-decay rates for a neu- 
tral atom placed in the vicinity of a dielectric half-space. 
We successfully re-derive the well-known results of phe- 
nomenological methods and broaden them by providing 
the asymptotic expansions that quantify the influence 
of absorption on the standard Casimir-Polder force cal- 
culated in Ref. [l^. We use only standard methods of 
quantum field theory, in a similar way as this is done in 
condensed matter theories. This requires the calculation 
of quantum propagators, most notably that of the elec- 
tromagnetic field. We show that this task is non-trivial 
but manageable. Inspired by the results of Ref. we 
find an exact solution of the Dyson equation satisfied by 
the photon propagator. In Appendix [C] we make con- 
tact with the phenomenological noise-current approach 
and calculate the photon propagator using the electro- 
magnetic field operators constructed on the basis of the 
noise-current operators Q. 



II. CONSTRUCTION OF THE MODEL AND 
HAMILTONIANS 

We are aiming to study the electromagnetic interaction 
between a quantum system, e.g. an atom, and a macro- 
scopic absorbing dielectric body. To this end we use the 
model of absorbing dielectrics developed in Q but gen- 
eralized to inhomogeneous dielectrics. The dielectric is 
modelled by a continuum of quantized harmonic oscilla- 
tors " the polarization field. This, in turn, is coupled 
to another set of quantized harmonic oscillators - the 
reservoir, the presence of which leads to damping in the 
polarization field so as to allow the absorption of radi- 
ation. These coupled quantum fields interact with the 
electromagnetic field via the coupling of the polarization 
field to the electric field. It turns out that the subsystem 
consisting of the reservoir, the polarization and the elec- 
tromagnetic field is exactly soluble, at least for simple 
geometries of the dielectric. Therefore, the interaction of 
the atom with the dielectric can dc facto be reduced to 
the interaction of the atomic dipole with the 'dressed' 
electromagnetic field, that is the electromagnetic field 
corrected for the presence of an absorptive body. This 
approach builds on the theory developed in [2l| where 
the interaction between an atom and a point-like absorp- 
tive dielectric (i.e. damped harmonic oscillator) was ad- 
dressed. The crucial difference is that for description of 
the interaction with a point-like absorber the 'dressed' 
electromagnetic field is required only perturbatively, but 
in the case of an extended absorbing body one needs 
to find the 'dressed' electromagnetic field exactly if one 
wants to capture accurately the interaction with an atom 
or other quantum system. 

Our starting point is the Lagrangian density describing 
the complete dynamics of the electromagnetic field and 
the dielectric, 

^0 = -CeM + Cp + Cr + Cp-EM ■ (1) 

The various constituent parts are: 

(i) The Lagrangian density Cem of the free electromag- 
netic field: 

/:EM = |E2(r)-^B2(r), (2) 

where E(r)_is the electric field and B(r) is the magnetic 
induction |22| . 

(ii) The Lagrangian density Cp of the polarization field: 

Cp = ixX(r) - ^Mw^^X^ir). (3) 

The field X is the dipole moment density of the con- 
tinuum of harmonic oscillators describing the dielectric. 
The strength of the restoring force acting on the po- 
larization oscillators is determined by the combination 
A4u!^. Hence, for a fixed absorption frequency wt of 
the dielectric, the 'mass' Ai is the parameter that de- 
termines the susceptibility of the polarization oscillator 
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to an external agent. It has dimensions of (mass) 



(length)" 



(dipole moment density) . In fact, the 



quantity (AieQUj^)^^ will turn out to be the polarizabil- 
ity of the dielectric at zero frequency [l2j. The absence 
of derivatives with respect to r in Eq. ([3]) implies that 
the polarization oscillators at different points in space 
arc mutually independent resulting in a model with no 
spatial dispersion. 

(iii) The Lagrangian density £r of the reservoir, in- 
cluding its coupling to the polarization field: 



dv { ip.Y.(r) - ]^p,v^ [Y,(r) - X(r)] 



(4) 



The set of fields represent the dipole moment density 
of the bath oscillators at all bath frequencies v, and the 
parameter has dimensions of (mass) x (length)"^ x 
(dipole moment density)"^ x (frequency)"^. The cou- 
pling of the bath to the polarization field leads to the ap- 
pearance of term proportional to X(r, t) in the equations 
of motion for the polarization field; hence it is responsible 
for damping [2^, [23l (cf. also Appendix [B]). The 'masses' 
of the bath oscillators pi, vary continuously with index v 
and describe the strength of the coupling between a sin- 
gle polarization oscillator and the continuum of reservoir 
oscillators for different frequencies v. The precise profile 
of pn is chosen so that the desired absorption spectrum 
is obtained f25| . 

(iv) The Lagrangian density £p_em describing the in- 
teraction of the polarization field with the electromagnetic 
field: 



Cp-EM = -.9(r)X(r) • E(r). 



(5) 



The dimensionless coupling function g(r) specifies 
where the interaction is taking place, i.e. 



5(r) 



1 inside the dielectric 
outside the dielectric 



(6) 



Thus g{r) describes the geometric shape of the dielectric 
body and limits the interaction to its interior. There- 
fore it is inconsequential whether the polarization field 
X(r) is defined in the whole of space or only in the in- 
terior, but the latter would cause unnecessary technical 
complications later on. 

It is straightforward to identify the canonical momenta 
and obtain the corresponding Hamiltonian densities 



"Hem 
Hp 

Hp-K 
'Hp -EM 



-^D2(r) + -^B2(r), 



P2(r) 1 



2M 

QO 

dv 





7Wa;|X2(r), 



2 

Zg(r) 
2p. 



^P.-^Y^(r) 



poo 

/ di/p.z/2x(r) .Y,(r), 
Jo 

.9(1"), 



-D(r) • X(r). 



(7) 

(8) 
(9) 

(10) 
(11) 



For convenience we have separated out the polarization- 
field reservoir coupling Hp-R and also the part of the 
Hamiltonian that just shifts the eigenfrequency of the 
polarization field. 



Hs = ^ rd-P.-^X^(r) + J^X2(r). 



(12) 



The first term of ([T^ arises due the coupling between 
the polarization field and the reservoir, whereas the sec- 
ond term is caused by the coupling between the elec- 
tromagnetic and polarization fields. Equations ([71)~ (fT^ . 
accompanied by the set of the equal-time commutation 
relations 



[A(r),i?,(r')] 
[^.(r),P,(r')] 
[n.(r),%,,(r')] 



z;ie'J™V;,(5(3)(r-r'), (13) 
zM,,,5(3)(r-r'), (14) 



allow one to derive the equations of motion for the inho- 
mogeneous damped-polariton model, cf. Q. The dielec- 
tric displacement D(r) = eoE(r) -I- (7(r)X(r) is the nega- 
tive of the momentum conjugate to electromagnetic vec- 
tor potential A(r), as it should be [2^. This is assured 
by the correct choice of coupling ([5]). As already men- 
tioned, the Hamiltonian density T-Ls, Eq. (fT2)) . shifts the 
eigenfrequency wt of the polarization field, i.e. 



,~,2 



1 

M 



— / dvpyV 



eoM ■ 



(16) 



The second term contains the parameter that pertains 
to the shape of the absorption spectrum. For our choice 
of i/-depcndence (see Appendix [B)) . it turns out to be in- 
finite. However, this is not problematic as the equations 
of motion for the fields and hence all observable quanti- 
ties, most notably the dielectric function, stay finite and 
physically meaningful. Furthermore, the last term of Eq. 
(|16p in principle introduces a position-dependence of the 
frequency ujt through the coupling function g(r). While 
not yet apparent at this stage, this position-dependence 
will turn out irrelevant. Hence, for now we set g{r) = 1 
in the expression for the frequency shift but shall explain 
later on why we are allowed to do so. With that we can 
incorporate equation (jl2[) into the Hamiltonian density 
of the polarization field and write 



-Hp = 



P2(r) 1 
2M 2 



(r), 



with 



1 

M 



dvp^v^ 



(17) 



(18) 



where by hindsight, we have introduced the symbol Wp = 
^ in analogy to the plasma frequency in metals 



Our aim is to investigate the influence of an absorbing 
dielectric on the properties of an atom, such as the shifts 
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in its energy levels and spontaneous decay rates. The 
complete Hamiltonian, including the atom, reads 



H 



n 



EM 



-n 



P-R 



n 



P-EM - 



Ha- 



EM) 



(19) 



where "Ha is the Hamiltonian density of the atom and 
Ha-em describes its coupling to the electromagnetic 
field. We consider a one-electron atom and treat the 
atomic electron non-relativistically by representing it as 
a quantum of the Schrodinger field satisfying fermionic 
anticommutation relations. The Hamiltonian density Ha 
of the non- interacting atomic electron (i.e. in the absence 
of interactions with the quantized electromagnetic field) 
can be written as 



Ha = *^(r) 



-|lv2 + y(|r-R|) 
2 m 



*(r), (20) 



where \I'(r) is the Schrodinger field operator satisfying 
the anticommutation relation 



{*(r),*t(i.)} = ^(3)(i._. 



(21) 



and V{\r — R|) is the potential due to the immobile nu- 
cleus which we choose to be located well outside the di- 
electric (i.e. at least few Bohr radii away) at a position 
R, so that there is no wave-function overlap between the 
atom and the solid. The atom is coupled to the electric 
field via its electric dipole moment. The Hamiltonian de- 
scribing this coupling in the dipole approximation may 
be written as 



Ha^ 



EM 



-/x-E(R), 



(22) 



i.e. we evaluate the electric field at the position of the 
nucleus. Here /x is the electric-dipole moment operator 
which depends on the second-quantized fields and 
It is convenient to expand the field operator ^'(r) in terms 
of a complete set of atomic wave-functions satisfying 



-V2-hy(|r-R| 



(/)„(r) = S„0„(r). (23) 



2m 
If we write 

vl'(r)=^c™^™(r), M't(r)=5^cL</';.(r) (24) 

m m 

then it follows from Eq. ([?!]) that the operators c,„ and 
cj( satisfy the equal-time anticommutation relation 



(25) 



We use equations ([M]) and ([^5]) to rewrite the Hamilto- 
nians Ha and Hem in a more useful form, 



Ha 

Ha-em 



n 

-e^c|cj(i|p|i) 



E(R) 



(26) 
(27) 



where —e{i\p\j) are the dipole matrix elements. 

We shall follow the field theoretical approach of [2l| 
to calculate the energy-level shifts and modified sponta- 
neous decay rates. In order to do so we need to locate the 
poles of the atomic propagator, which, once interactions 
have been switched on, are accessible only pcrturbatively. 
For these perturbative calculations we need to work in 
the interaction picture where the general expression for 
the perturbative expansion of a Green's function of the 
field under the influence of the interaction Hi is [2^1 



n=0 



n+l 



dti 



dtn 



X {n\ T [vl/(r, t)*t (r', t')Hi{h) . . . i?i(i„)] .(28) 

^I* is now the field operator in the Heisenberg picture and 
|r2) is the exact ground state of the non- interacting sys- 
tem. The subscript 'conn' indicates that only connected 
diagrams contribute, as disconnected diagrams drop out 
in the normalization of |f2). 

Wick's theorem states that the terms appearing in the 
expansion (|28p . when written out explicitly for a specific 
interaction Hamiltonian, turn out to be given entirely in 
terms of the propagators of the non-interacting fields. We 
shall proceed by determining the non-interacting propa- 
gators of the atom, the polarization field, the bath, and 
the electromagnetic field. Then the interaction of the po- 
larization field with the reservoir shall be treated exactly 
to all orders. Once this is accomplished the correction to 
the electromagnetic-field propagator caused by the pres- 
ence of the absorptive dielectric can be calculated, which 
shall also be done exactly to all orders. This is going to 
give the dressed photon propagator that enters the final 
perturbative expansion of the atomic propagator whose 
poles give the energy-level shifts and changes in the tran- 
sition rates. 



III. UNPERTURBED FEYNMAN 
PROPAGATORS 

A. Atomic-electron propagator 

The unperturbed atomic-electron propagator corre- 
sponding to the Hamiltonian (|20p or equivalently ([26]) 
is defined as the time-ordered expectation value 

G(°)(r,r',i,t') = T [*(r, i)*t(i.', f)] (29) 

where 5' is the Schrodinger field operator in the Heisen- 
berg picture and jO) is the exact ground state of the 
non-interacting system. We substitute the field operators 
written in terms of the atomic eigenfunctions, Eq. (|24p . 
while remembering that we are in the Heisenberg picture 
where the operators c/ and are time-dependent, and 
find 



G(°)(r,r',i,<') = E'^'W 



(30) 
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with 



Giy,t')^^^{^}\T[ci{t)cl{t')]\n). (31) 

The time-dependence of the fermionic annihilation and 
creation operators is governed by the Hamiltonian (P^ . 

cm{t) - c,„(0)e-*^'"*/^ 4(i) = 4(0)e*^'"*/\ (32) 

With that we can determine G|"'(i,t') and obtain 

GT2it t') = - - i')e-'^'(*-*'^/"<5,„., (33) 

where we have used the definition of the time-ordering 
operator and the fact that the vacuum state is an- 
nihilated by Cm(0). Since G[2it,t') is in fact dependent 
only on the time difference t — t' , we can work with its 
Fourier transform with respect to t — t' = t 



Gfl{E)=. / dre-^/^Gj^(r) = 



1 



E - El +17] 



(34) 

With this convention of Fourier transformation the irj- 
prescription ensures the correct causal behaviour of the 
propagator and guarantees the convergence of the inte- 
grals. 



B. Photon propagator 

To calculate the zcroth-order propagator of the dis- 
placement field D(r,t) whose dynamics is governed by 
the Hamiltonian ([7]), which we emphasize does not in- 
clude the coupling term (fTTj) . we note that the Heisenberg 
equations of motion imply 



dD,{v,t) _ 1 



dt fiQ 



gl £ Vji;fe(r,i), 



(35) 
(36) 



where e'-''^ is the Levi-Civita symbol and the sum over 
doubly occurring Cartesian indices is implied. Thus the 
displacement field D(r, t) satisfies the homogeneous wave 
equation 

(V,V, -%V2)i5j(r,i) +Aioeo^A(r,<) =0. (37) 
The formal definition of the photon propagator reads: 

A,(r,r',t,i') = -^(0| T[A(r,i)^j(r',0] |0), (38) 

where Di{r,t) is the displacement field operator in the 
Heisenberg picture and |0) is the exact ground state of 
the non-interacting electromagnetic field. We proceed by 
applying the differential wave-operator that appears in 
(|57|l to this definition of the propagator, but we need to 



take care when applying the time derivative to a time- 
ordered product and observe that 



^T[A{t)Bit')]=Sit-t')[A{t),Bit)]+T 



dAjt) 
dt 



Bit') 



Thus we find that the displacement field propagator 
Df^\r,r' ,t^t') satisfies the following differential equa- 
tion: 



A-t' 



(27r) 



'-^S{t - t') jd\{mk - <5,,.q2)e^i ("--'), (39) 



where we have used the commutator 



dt 



:Dk{r\t) 



-(V,V.-5.fcV2)<5(3)(r-r') 
Mo 



and the fact that spatial derivatives commute with time- 
ordering operator. From Eq. ([39|) it is clear that the free- 
space photon propagator is translation-invariant in space 
and time, i.e. it depends only on the differences r — r' 
and T = t — t'. Therefore one can find the solution of 
the differential equation through Fourier transformation. 
First wc note that Maxwell's equation ([55]) implies that 
the displacement field is transverse, so that its propagator 
satisfies 

V,Df,\r-v',t-t')^0. (40) 
Introducing the Fourier transform of the propagator 

dII^ (q, Lo)^ jcP{T~ r')e-*'i-('-'"') 

/•OO 

X / d(^-i')e*"^*"*'^-Dlfc'(r-r',t-t') (41) 



(42) 



we readily obtain its spectral representation 
(0)/ . . .^ , kkG^ - (Mk 



^2 _ q2 _^ iry 



We have displaced the poles in the denominator by ir] 
so that D^^i}{r — r',t — t') has the appropriate causality 
properties of a Fcynman propagator. 



C. Polarization field propagator 

The Hamiltonian density (fT7|l describes a collection 
of mutually independent harmonic oscillators. The fact 
that the harmonic oscillator at r is unaffected by the 
oscillator at r + dr allows us to introduce creation and 
annihilation operators, bt(r) and b(r), for each harmonic 
oscillator, 



X(r) 



2A^wt 



[bt(r)+b(r)] 



P(r)=n/^ [bt(r)-b(r)] 



(43) 
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The operators b'''(r) and b(r) satisfy the equal-time com- 
mutation relations 



6,(r),6](r)l =5,,<5(3)(r-r') 



(44) 



which follow directly from their definition and Eq. p4p . 
The operator 6i(r) annihilates the ground state of the 
oscillation in the i-th direction at r. Using this property 
together with the commutation relation (|44p . we can di- 
rectly evaluate the polarization field propagator defined 
as 



--{n\T[X,{r,t)X,{r',t')]\n}. (45) 



Here Xi{r,t) is the polarization field operator in the 
Heisenberg picture and |f2) is the exact ground state 
of the non-interacting polarization field. When written 
in terms of the creation and annihilation operators, the 
Hamiltonian density (fT7|) is of course diagonal in b^'(r) 
and b(r) so that the time dependence of the creation and 
annihilation operators is harmonic, 

b(r,i) = b(r,0)e-''^'^*, b^(r, t) = b1'(r, 0)e^'^'^*. (46) 

We substitute the polarization field operators ex- 
pressed in terms of the ladder operators into Eq. ((45l) 
and observe that, due to the orthogonality of states, only 
terms proportional to bibl contribute. Taking care of 
the appropriate time ordering of operators and using the 
commutator p4| to move any annihilation operators to 
the right of creation operators, so that they act on the 
vacuum state we readily obtain 

X \e{t - t')e-»'^T(t-t') ^ 0(^f' _ ^^^+^u■r{t-t')^ ^ (47^ 



with the frequency Cjt as defined in Eq. (|T8|l . We shall 
need the Fourier transform of the polarization propagator 
with respect to the time difference t — t' 

/oo 
d(i-t')e*"^*"*'^<°^(r-r',t-t'), 
-OO 

(48) 

which is easily obtained from Eq. (I47p and reads 



1 



— %<5(3)(r-r'). (49) 



Since the polarization field operators are mutually inde- 
pendent there is no need for any special consideration 
of the boundaries of the dielectric medium at this stage. 
The boundaries are taken into account through the cou- 
pling function g(r), Eq. ([S]), and the field equation for 
the electromagnetic field, once coupled to the polariza- 
tion field, will include the physical processes of refiection 
and refraction as it should. It is worth pointing out that 
an artificial restriction of the free polarization field to just 
the interior of the dielectric would lead to a much more 
complicated free propagator thereby causing unnecessary 
technical complications while not describing any different 
physics. 



D. Reservoir propagator 

The dynamics of the non-interacting reservoir field is 
governed by the Hamiltonian © which describes a set 
of independent harmonic oscillators. The propagator for 
the free reservoir field can be obtained by repeating the 
same steps as for the derivation of the propagator for 
the free polarization field in Section IHI CI Therefore, we 
do not repeat the details of the derivation but simply 
point out the similarity of the structure of the result to 
Eqs. (|47p and P^)) . In the time domain the reservoir 
propagator reads 



H^°\r -r',t- t', V, v') = --^^.^.^(s) (r - r')6(v - v') 
X \B{t - tOe"""^*^*'^ + ^(^' - <)e+"'**"*'^l • (50) 



Its Fourier transform with respect to t — t' is given by 



d{t ~ t') e'"(*-*')i?^)(r - r', t - t', i., z/'), 

o 



(51) 



IV. DRESSED PROPAGATORS 



Having gathered all the unperturbed propagators, we 
can proceed to work out the propagators for the coupled 
fields. We are going to use a diagrammatic approach to 
illustrate the workings of perturbation theory, i.e. we 
represent each term of the perturbative expansion (|28p 
by an appropriate Feynman diagram (e.f. e.g. psj). To 
proceed with that, we need to lay down the Feynman 
rules for our approach. We have four different free prop- 
agators; accordingly, we associate with them four distinct 
lines: 



= ihGf{t,t') 



II 



r', t' ran r, t 



We shall need to consider three interaction Hamiltonians 




7 



in turn, i?p-R, -ffp-EM and Ha-em- 

Hp-EM = -- /d3rg(r)D(r) • X(r), 
fo J 

Ha~em = V c\cj ■ D(R). 

^0 ~r 



(52) 
(53) 
(54) 



Note that we have introduced the shorthand fi^j = 
{i\fi\j) for the matrix elements of the atomic electric 
dipole moment operator fi. These interaction Hamilto- 
nians yield the following Fcynman rules for the vertices 
between the lines defined above: 



kl 



33 



kl 



r. t 



ri. ti 



h>,pq 



To compute a diagram one has to sum over all inter- 
nal indices and integrate over internal times, internal co- 
ordinates, and reservoir oscillator frequencies u and ly' . 
As mentioned earlier, the subscript 'conn' in Eq. ([28]) 
means that the summation in that equation runs only 
over those terms that correspond to connected Feyn- 
man diagrams. Furthermore, topologically equivalent di- 
agrams, i.e. those that can be obtained by permuting 
the factors Hi{ti) in Eq. p8|) . are counted only once, and 
therefore we have omitted the factor of 1/n! that would 
otherwise have arisen in a straightforward expansion of 
the time-ordered exponential of the interaction Hamilto- 
nian in perturbation theory. 



A. Dressing the polarization line 

The polarization field interacts with the reservoir; all 
these interactions in their entirety "dress" the polariza- 
tion field. We choose to represent the dressed polariza- 
tion propagator by a bold dashed line: 



= ihKmJr.v'AA') 



mn 



From the interaction Hamiltonian (|52p and the associ- 
ated the Feynman rules one can see that the polarization 
line can only ever connect to exactly one reservoir line. 



Hence the complete set of all possible interactions cor- 
responding to the expansion ([28]) is represented by the 
following sequence of Feynman diagrams: 

+ • • 

+ • • • • 

+ ... 



The equivalent analytical expression is the Dyson equa- 
tion for the dressed polarization propagator; it reads 

i^„„,(r,r',i,i') = Al°^(r,r',t,t') 

/OO 
i,p 

X kI^} (r, ri , t, t^)H[l^ (ri , , ti , ^2 , i^, v')Kpn {^2 MX) 

Despite being an integral equation, the above equation 
is easily solved exactly. Upon substituting Eqs. ((49l) 
and ()5ip we can easily carry out the spatial integrations. 
Then we Fourier transform with respect to t — t' , 



dt2 l(Pvi l(Pr2 / dv I dv' 

"'0 



r,r ]uj] 



d{t-t')e''^^*-''^K„,n{r,r',t-t') , 



and find the following expression for the dressed polar- 
ization field propagator in the frequency domain 



(55) 



with 



M Jq '^^ U;2 - iy2 ^ jy^ 



(56) 

Note that K{uj) is an even function of u. The plasma 
frequency wp was defined below Eq. (fT8)) . 



B. Dressing the photon line 

The coupling ([55]) between the dressed polarization 
field and the electromagnetic field has formally the same 
form as the coupling ([52]) between the bare polarization 
field and the reservoir. By analogy with the previous 
section, we write down the graphical equation for the 
dressed photon propagator as 



*\/\/W\, 



where the bold wavy line denotes the dressed photon 
propagator i.e. 

^\/\/\/\/\/\/\, = ifiDki{r, r', i, t') 
r', t' kl r, t 
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The corresponding analytical expression reads 
D,k{r,r',t,t') = D^l\r-r',t-t') 

+ ^^1/ dt2 d^Yi d^Y2g{Ti)g{Y2) 

Jt'' 

X (r - r 1 , t - 1 1 ) A^,7 (r 1 - r2 , 1 1 , t2 ) Afe (r2 , r' , t2 , t') ■ 

(57) 

Now recall the discussion following Eq. (fT6|) of the shifted 
eigenfrequency a)x of the polarization field. It enters the 
Dyson equation (|57|) through the dressed polarization 
field propagator Kji{Yi — Y2,ti,t2)- As we noted ear- 
lier, according to Eq. ([T6)) . the shifted eigenfrequency wt 
suddenly jumps at the boundary of the region where the 
polarization field interacts with the electromagnetic field 
i.e. where the coupling function gij) = 1. However, it 
is now apparent that this discontinuity is unproblematic 
because all spatial integrations in Eq. ([57]) are limited to 
the region of space where g(r) = 1. 

To simplify Eq. ([57]) we note that one of the spatial 
integrations is trivial due to the 6 function in the dressed 
polarization field propagator (|55p . Fourier transforming 
with respect to the time difference t — t', 



direct iteration. The iteration method that we shall em- 
ploy is inspired by Ref. In order to explain it, we 
write Eq. ([5^ symbolically as 

D = D° + KD° ® D. (61) 

Iteration of this equation yields the expansion 

D = D° + KD" + K^D^ ® D° ® D° + . . . , (62) 

which proves especially useful if the action of the operator 
O = on the free-space propagator amounts to a 
simple multiplication, i.e. if 

OD° ^ D° ® D" ^ CD°, (63) 

where C is some constant. Then Eq. (f62|) becomes a 
geometrical series 



D = D" {1 + KC + K^C^ + K^C^ 



(64) 



which wc know how to sum to all orders. 

An alternative approach, which we sketch in Appendix 
lAl consists of converting the integral equation ([59|) to a 
differential equation supplemented by Maxwell boundary 
conditions. In addition, in Appendix [Cl for comparison 
with other theories, we construct the photon propagator 
using yet another, completely different method based on 
the phenomenological noise-current approach of Ref. [1] . 



D^k{Y, 



)= / d(t-i')e^"'*"*'^Afc(r,r',t,t') , (58) 



we find the Dyson equation for the dressed photon prop- 
agator in the frequency domain: 

Afe(r,r';w)=i?(°)(r-r';c.) 



Jd^YigiYi)D[fiY - Yi;uj)D,k{YuY';u). (59) 



Here -R'(aj) is a complex- valued function of frequency that 
has originated from the dressed polarization field propa- 
gator and is given in ([55)1 : it will be shown to be related 
to the dielectric permittivity. Note that the dimension- 
less coupling function g{Y) describing the geometry of 
the dielectric medium, as defined in Eq. ([6]), is the only 
way the geometry enters in the calculation, by effectively 
defining the limits of the spatial integration in Eq. (|59p . 

^ik (i" ~ is the free-space photon propagator in co- 
ordinate representation i.e. the inverse Fourier transform 
of Eq. dUl), 



dII\y^y\c.) 



(27r)3 



iq(r-r') ^»fcq^ - gigfc 



IT] 



(60) 

The solution of the integral equation ([59]) is much less 
trivial than that of the equivalent equation for dress- 
ing the polarization line in Section llV Al This is because 
translation invariance is lost when an inhomogeneous di- 
electric is introduced into the system. Here we report 
two ways of tackling the problem. First, we demonstrate 
that it is possible to solve the integral equation ([5^ by 



z 







FIG. 1. The atomic dipole is located a distance Z away from 
the dielectric half-space of complex and frequency-dependent 
permittivity e(tj). The transverse propagator Difc(r,r';cj) of 
the dielectric displacement field in this geometry is given by 
Eq. (f96l) . 



Let us now concentrate on the example geometry of 
a dielectric half-space occupying the z < region of 
space, cf. Fig. [U for which the coupling function g{Y) 
in Eq. (|59|) becomes g{Y) ~ 9{~z). Due to the bound- 
ary the problem has lost translation invariance in the 
z-direction, but not in directions parallel to the surface. 
In other words, the propagator depends only on the dif- 
ference r|| — rj|, but separately on z and z'. It is con- 
venient to work with quantities that have been Fourier 
transformed with respect to ry — ; e.g. for the dressed 
photon propagator we have 

D,,iz,z') = |rf2(r,| -r||)e-'iH--„-r;,)^^^.(r|| _r||,z,z'), 

(65) 
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where for notational convenience we have suppressed the 
dependence on q|| and w. Once Fourier transformed with 
respect to r|| — rj|, the integral equation (|59p becomes 



D,k{z,z') = D'i;j{z-z') 



dz^DfMz- zi)D,k{zi,z').{m) 



As is easily checked, this does not lend itself to iteration 
as it stands. Following Ref. [1^ we introduce an addi- 
tional integral equation in order to enable the iteration 
process: 



D,k{z,z') = Df^{z-z') 



K{u) 



/ dziD^f {z- z,)Djk{zuz'). {67) 
Jo 



Here D^^ {z — z') is the Fourier-transformed photon prop- 
agator in a bulk medium, i.e. the solution of Eq. ([5^ 
with g(r) = 1. In order to justify Eq. (jHT)) let us recall 
that the part of the Hamiltonian density that describes 
the interaction of the photon field with the polarization 
field has the form 

H=Ho-^^X(r).D(r), (68) 

where 'Ho is the Hamiltonian density of the non- 
interacting electromagnetic field. Using the fact that 
0{—z) + d{z) = 1 we can also write 



1~L — l~Le 



9{z) 

(0 



X(r)-D(r), 



(69) 



where He ^T-Lq — X(r) • D(r)/eo is the Hamiltonian den- 
sity of the electromagnetic field interacting with an un- 
bounded dielectric. Therefore, we have a choice: we can 
either correct the free-space photon propagator for the 
presence of the dielectric half-space or, equivalently, cor- 
rect the photon propagator in a bulk dielectric for the 
absence of the medium in the other half-space. In other 
words, it is entirely up to us which Hamiltonian we take 
as the zeroth-order (exactly solvable) Hamiltonian when 
applying perturbation theory. The integral equation (j67p 
corresponds to treating the electromagnetic field interact- 
ing with bulk medium as the zeroth-order, solved part of 
the problem. 

To proceed we need to find D^^\z~ z') and Z?,-^' {z — z') 
appearing in Eq. ((66l) and ((67|) . which can, in fact, be 
read off from the representations of these propagators 
as two-dimensional integrals over the momenta parallel 
to the surface. To find Dfl}{z — z') we carry out the 
qz integral in Eq. ([60|) using the residue theorem. The 
result, written in a compact form, is 



^!°'(r-r',a.) 



2(27r)^ 



(V.Vfe - ^.feV^ 



d2q||g»qr(r||-r;,) 



Jk^\z-z'\ 



(70) 



where kz is the z-component of the wave vector in vac- 
uum and is given by kz = ^a;^ — q| + The square 

root is taken such that the imaginary part of kz is always 
positive. Eq. ([70]) shows that some components of the 
Fourier transform of the free-space photon propagator 
are singular when crossing the z = z' plane. 

For deriving the photon propagator D^'j}{r — r';uj) in 
a bulk medium, we set g{r) = 1 in Eq. ([55]) and then 
Fourier transform this equation with respect to r — r', 



IT] 



K{uj) 



eo 



ir] 



^r^q,^)- (71) 



This matrix equation becomes an algebraic one when one 
takes the transversality of the propagator, qjO')^^ (q, w) = 
0, into account. The calculation is straightforward and 
in coordinate space we obtain 



eo^(^) 
(27r)3 



(72) 



Note that the function ^(w) that appears in Eq. (jTS]) 
should not be interpreted as the dielectric function of 
the bulk medium. It is an even function of the frequency 
cj and may be written explicitly as 



e(^) 



eoTW 



2 2 ^ / , 
M Jo 1 



(73) 

-1 



where we have used Eq. (j56p . Thus the w-dependence 
of ^(w) is not consistent with the causality requirements 
usually imposed on response functions, i.e. with Kramers- 
Kronig relations. This is because we have calculated a 
Fcynnian propagator and not a retarded Green's func- 
tion. The dielectric function of this model is discussed in 
Appendix |b1 We derive the required Fourier representa- 
tion of the propagator in the bulk medium by carrying 
out the qz integral in Eq. ([72)1 and obtain 



^eo^(^) 
' 2(27r)2 



X d qye 



'2n„„»q||-(r||-r;|) 



aik-d\z-z'\ 



(74) 



in complete analogy with the formula for the free-space 
propagator, Eq. ([70]) . Here kzd 



/^(aj)w2 — q| is the 

z-componcnt of the complex wave vector in the medium 
with an always positive imaginary part. 

We may now proceed by substituting Eq. (|66p into Eq. 
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mi 



'^^ rdz,Dlfiz-z^)D^iz,^z') 
Jo 

^ Hdz, f dz^D^fiz-z,) 

Jo J-oo 



xDfi\zi-Z2)Dikiz2,z') (75) 

and focusing our attention on the solution of the case 
z < and z' > 0, i.e. when the source is located out- 
side the dielectric and the observation point is inside the 
material. The solution for the case z,z' > can then 
be obtained by applying the integral equation (|66|). The 
advantage of introducing Eq. ([75]) is that it facilitates it- 
eration as it turns out that when Dik on the right-hand 
side is replaced by D^^ the action of the double-integral 
operator in the last term reduces to a matrix multiplica- 
tion, 

roo 1*0 

/ dzi dz2Dlfiz-zi)Dfl\z,^Z2)DllHz2-z') 

Jo J ~oo 

^C,,D\k{z-z'),{n) 

with the matrix Cij independent of z and z' . In order 
to efficiently verify and make use of assertion (|76p let us 
point out some useful facts. First we recall that 

qH,u = [eP(q)e™(q) + e™(q)e™(q)] , 

(77) 

where q = (q|| , fc^) is the wave vector in vacuum and we 
have introduced the polarization vectors 



1 



e (q||)=| — r(-'7j/,fe,0), 

Iqiil 



e™(q||,M = 



(qii,fc^ 



Iqiil^ 



{Qxkz,(lykz, qn 



(78) 



1 



-{qxkzd,qykzd,-c^\) 



|qillV?Hw 

We have listed e^^{<\\\,kzd) explicitly to point out the 
additional factor of ^"^/^(w) in its normalization. In the 
following we will suppress the insignificant dependence of 
the polarization vectors on q|| . Relation ()77p is simply a 
statement of the completeness property of the polariza- 
tion vectors (l78l). but it allows us to write 



Lo e 



X 



z > z' 

eK-kz)el{-kz)e-'^'^'-'-'\ z < z' 



We emphasize that the above representation of the free- 
space propagator is not valid at the point z = z', where 
the z-derivatives in Eq. (f70|) acting on e^^^^l^^^ I would 
produce additional terms proportional to a delta func- 
tion. Similarly we have 



IK. 



A 



e\(kzdyl{kzd)e'^^-^'-''\ z>z\ 
eH-kzd)eU-kzd)e-''''-^''''\ z < z' . ^ 



Eqs. (|80p and (f8T|) show that the free-space and bulk- 
medium propagators can be split into separate contribu- 
tions from the transverse electric and transverse magnetic 
polarizations 



(82) 



Most of the further calculations are very much simplified 
if one takes into account that scalar products of polar- 
ization vectors with different z-componcnts are diagonal 
in the polarization indices, i.e. we have 



(83) 



The function / is equal to 1 for the TE mode, and for 
the TM mode it reads 



f'''\qz,Pz) 



QzPz 



(lft+Qz\hft +pI 



(84) 



This is very useful because it shows that not only a single 
propagator, as in Eq. ()82p . but also a product of propaga- 
tors can always be split into separate contributions from 
the transverse electric and transverse magnetic modes, 
i.e. we can always write 

. . .Dij{z - Zl)Dj;(zi - Z2)Dlk{z2 - z') .. . 

= ^ . ..Dx.ijiz - zi)Dx,ji{zi - Z2)Dxm{z2 - z'). . . . 

A 

This is true for an arbitrary number of propagators. 

We can now proceed to verifying Eq. ([75)1 . First we 
note that the arguments of all three propagators entering 
Eq. ([76|) have a definite sign. Indeed we have 



Zi < 0, Zi — Z2 > 0, Z2 — z' < 0. 



(85) 



(79jhus, it follows from Eqs. (pT)) and ([5T|) that the propa- 
gators entering the integral in Ec^. (|76p are given by 



so that the partial Fourier transform of the free-space 
propagator ((70|) may be written as 



A?(z - z') = _r^e-'^^-^^^ > ' et{kz)e^ikz), 



2k, 



= _!f^f,^k,{z-z')^pX( 
2kz ^ 

ieoC^(w)a;2 



-ik^d{z-z') 



2k, 



z > z' 



eHkz)e^ikz)e'''^^'-''\ 
e^{-kz)el{-kz)e-^^'^'--'''\ z < z' 



,(80) 



x^e^^(-Met(-M- (86) 

A 
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With this we can evaluate the integrals in Eq. (|76)) and 
find that 



dz, I dZ2Dlfiz - Zl)D'^;>{zi - Z2)D\1>{Z2 ~ z') 



(0), 



2 
A 



where we have used Eq. ((751) . Here rx is the Fresnel 
coefficient for reflection from a half-space. Since all the 
Fresnel coefficients for reflection and transmission at a 
half-space will be needed later on, we list them here: 



?'TE 



^TE 



kz + kzd 

2kz 



rTM 



tTM 



^{uj)kz + kzd 
2^/^kz 



i{Lo)kz 



The significance of Eq. ([87|) is that it allows us to iterate 
the integral equation ([75]) along the lines of Eqs. (|63l) and 
()64[1 . The iterative process is now straightforward and, 
thanks to relation ((87)) . leads to two separate geometric 
series for the two polarizations 



D\,ij{z,z') 



^ / dz,D[%iz^z,)Di%i^i-^^') 

JO 



1 - 



1-rl 



1-rl 



(89) 



These geometric series can easily be summed up to all 
orders to give the exact photon propagator for the case 
z < 0, z' > 0. In order to cast the result into a familiar 
form, we explicitly evaluate the integral in the second 
line, which requires some care. The integral that needs 
to be evaluated is 

z') = ^ rdz,D%{z - z,)d'-^],{z, - z'). 
^0 Jo 

(90) 

Here the argument of D'^^'^ is always negative, z — zi < 
0, whereas the sign of zi — z' can be both positive and 
negative. Therefore we need to take into account that the 
propagator (|70p is discontinuous at zi = z' and contains 
singular terms proportional io 5[zi — z'). 



In order to correctly evaluate the integral ((90)) wc rep- 
resent the differential operator in Eq. ()70|) using the po- 
larization vectors written out in terms of derivatives. Us- 
ing the completeness relation of the transverse polariza- 
tion vectors we may symbolically write 

V.V, - J.feV^ = -V^ ^ e^(V)e^(V). (91) 



With this, the propagators entering the integral ()90p are 
given by 

xgf^.ki-^'l^ (92) 



2k 



zd 



eti~kzd)en~kzd) 



xe 



(93) 



Note that in we have changed the z-derivatives to 
act on z' rather than on z so that they could be pulled 
outside the integral in ((90)) . Now it is straightforward to 
demonstrate that the integral ((90)) is given by 

I>^,{z,z') = D[%iz-z') 
2kzd t\ 

whose first term exactly cancels the bulk dielectric prop- 
agator in the first line of Eq. ((89)) . The remaining term 
yields the final result 



Aj(2, z') 



^2k7 



J2 [a^)efhkzd)e^i~kz)tx] 

A 

xe-''''-''+'''''' , (95) 



with the transmission coefficient as given in ((88)) . This 
formula describes the vacuum-dielectric transmission i.e. 
it is valid for z < 0, z' > Q. It is a straightforward cal- 
culation to plug equation ((95)) into ((66)) and obtain the 
photon propagator for the case z,z' > 0. In the region 
z' > the final result for the dressed photon propaga- 
tor Fourier transformed back to coordinate space may be 
written as: 



A, (r, r';u)= e{z)Dlf (r - r'; ..) - ^ ^ | d^q,, ^^e^''" -(r,, -1, ) 

x|0(-z) [^(w)e,^(q||,-Me,'(q||,-fc.) ^a] 6"^'=-^+'^'=^^' + 0(z) [e^(q|| , fc,)e^^(q|| , -fc.) rx] e^'=^(^+^')|. (96) 

I 



In the calculations of the energy-level shifts of an atom placed outside an absorbing dielectric material, to be dis- 
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cussed in the following section, we shall need the prop- 
agator for the case z, z' > 0. In that case Eq. (|96|) 
shows that the propagator splits into a free-space part 
-D-°^ (r— r'; cj), which is not interesting as it just yields the 
standard (position-independent) Lamb shift, and a cor- 
rection due to reflection at the boundary, which we shall 
call Dl'^\r,r';uj) and which gives a rise to the position- 
dependent Casimir-Polder shift. As we shall treat the 
atom-field interaction in the dipole approximation, we 
are going to need the reflected part of the propagator 
-D,[j'' (r, r'; w) evaluated at equal arguments r = r' = R, 
where R = (0, 0, Z) is the position of the atom. In that 
case it simplifies considerably and can be written in the 
form 



D('-'(2;w) 



8^ JO 






^11 F 







2(zf|r™ 



(97) 



with kz = Y ^ 'Zy + s-s before. Note that we 
have gone to polar coordinates, qx = q\\ cos (/>, qy = 
q'l|sin0, where the azimuthal integration annihilated 
the off-diagonal elements of equal-argument propagator 

As a final remark we would like to comment on the 
convergence of the series in Eq. ([SQ]) . It clearly converges 
provided 



1 



< 1. 



(98) 



However, there does not seem to be a physical signifi- 
cance to this condition. That the result for the propa- 
gator can be extended by analytic continuation to wave 
vectors not satisfying the condition ([98]) can be shown by 
solving the corresponding boundary-value problem, the 
procedure for which we sketch in Appendix |^ 



V. ATOMIC PROPAGATOR AND ELECTRON 
SELF-ENERGY 

In order to investigate the perturbative expansion of 
the atomic propagator ([28]) . we use the expansion in 
terms of atomic eigenstates, Eq. ((24|) . and then work 
with the atomic propagator in that basis. In analogy 
to Eqs. (pO|l and (|3ip for the unperturbed propagator, 
we obtain 



n=0 



n+1 



dti 



dtr, 



X ( Vt 



[c,;(t)c|(t')i^A-EM(il) . . . -ffA-EM(t„)l ^) (99) 
L J / conn 

By using Wick's theorem to evaluate the ground-state ex- 
pectation value of the time-ordered product of operators. 



one easily sees that the zeroth-order term is a propagator 
for non-interacting system and that the first-order correc- 
tion vanishes because it is not possible to contract all of 
the operators. Therefore the lowest-order non-vanishing 
perturbative contributions come from terms of order e^. 
Diagrams to this order have two vertices and therefore 
include a pair of disconnected tadpole diagrams, which 
arc irrelevant as they go away in the process of normal- 
ization, and the physically important self-energy diagram 
which contains all the information about the energy-level 
shifts and decay rates. 




For calculating the perturbatively corrected, dressed 
atomic propagator it is in fact convenient to perform a 
partial summation and consider the following series of 
diagrams 



where the thick solid line represents the dressed atomic 
propagator. The location of poles of the such constructed 
propagator is much more straightforward to work out 
than for the propagator with strictly only one-loop cor- 
rections. The above is a graphical representation of the 
Dyson equation, which expressed analytically reads 



G,,{t,t') = Gf{Ut') + 



ih 



^2 / J r^mir^tm 
k.l.m 



— oo J — oo 



dti / dt2 



xG[°^ {t, h)Gl±{h,t2)DkiCR, R, h,t2)Gu{t2, t'). (100) 

We note that here and in the following the indices k and 
I label just Cartesian components, but the sum over m 
is a sum over intermediate atomic eigenstates, and i is 
the atomic state whose energy shift we are seeking to 
determine. To this end we Fourier-transform (|100p with 
respect to t — t', along the line of Eq. ([34| and find the 
dressed atomic propagator 

/oo 
d(t_t')e»(*-t')£AG,,(t-f') 
-OO 

= P p ^ ^ ^^777 (101) 

with the self-energy insertion 



ih 

27rel — 

" k,l,m 



i5fez(R, R;w) 



E — Tiuj — Em + irj 
(102) 

The self-energy insertion (|102p contains the dressed pho- 
ton propagator which in the case of an atom outside a 
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dielectric half-space comprises: (i) the free-space pho- 
ton propagator D^^-* (R, R; w), which gives a rise to the 
position- independent Lamb shift, and (ii) the reflected 
part £'^j^(R, R; a;), which yields the position-dependent 
Casimir-Polder shift. Thus the shift in the atomic energy- 
levels, given by the poles of Eq. (|101[) . can be written 



(103) 



As we want to work out changes in the energy levels al- 
ready corrected for the coupling between the atom and 
the free-space electromagnetic fields, we renormalize the 
energy-level shift by subtracting the self-energy associ- 
ated with the free-space electromagnetic field S^^"* {£) and 
consider 



We use the symbol Ei = Ei + ) to represent the 

atomic energy levels already corrected for the free-space 
Lamb shift and decay rates. 

Prima facie it may seem difficult to extract the energy 
shift from Eq. (jl04[) because it is an implicit equation 
whose right-hand side also depends on £. However, if 
the energy shift we calculate is small compared to the 
difference in energy between the state under considera- 
tion and its nearest dipole-connected neighbours (which 
it needs to anyhow for perturbation theory to be appli- 
cable) , then the shift can be extracted from Eq. (|104p by 
a single iteration leading to 



(104) 



Zj[;^(R,R;c.) 

Wmi + CJ - ij] 



(105) 



where we have abbreviated ujmi ~ — ^i- The w- 
integral in (|105p can be restricted to the positive real 
axis by writing 

(106) 



and noting that £'|.'.''(R, R; w) is even in w (see Section 
IIVBI and Appendix |B]) . Then the term proportional to uj 
is odd and vanishes when integrated over the real w axis. 
As the photon propagator is analytic in the first quad- 
rant of the complex cj-plane, it is permissible, provided 
> 0, to rotate the contour of w-integration by 7r/2 
i.e. w — >■ icj. This applies when one considers an atom in 
its ground state. 

However, for an excited state i of the atom one has 
'^mi < 0, which means that there will be poles in the 
first quadrant of the w plane due to the denominator in 
(|105p . We would also to remark that the Fresnel reflec- 
tion coefficients have poles in the complex plane at the 
location of trapped electromagnetic modes, which is not 
an issue in the case of a dielectric half-space but arises 



e.g. for a dielectric slab [29| and other systems capable 
of wave- guiding [30 |. 

We recah that £'^^^(R, R; cj) is diagonal, cf. Eq. (P7)) . 
and write down the final result for the energy shift in the 
form 



(107) 



with A£i and A5* given by 



,k |2 



duj- 



2 , ,,,2 ^kk 



i?^y(R,R; iw) 



(108) 



= i E lMLl'^i;^(R,R; |w™|)^(-^™) (109) 



k,rn 



where = |(m|/z'^|z)| are the matrix elements of the 

A:-th component of the electric dipole moment operator. 
The quantity A£* is the contribution to the self-energy 
that originates from the poles in Eq. (jl05|) that arise for 
an excited state i for which ujmi < 0. Expressions equiva- 
lent to Eqs. pOSp and p09|) have been derived before by 
different methods, e.g. by linear response theory [sil . l32j 
or later by the noise-current approach to phenomenolog- 
ical QED H. 

The shift A£i is real because it is a convolution of 
atom and field susceptibilities which are real at complex 
frequencies [3l| . However, AS* is complex; its imagi- 
nary part modifies the decay rates of excited states. In 
summary, we have 



AEi = A£, +Re(Af*) 
AF,; ^ -|lm(A£*) 



(110) 



where AEi are the rcnormalized energy-level shifts and 
ATi are the changes in decay rates. 



VI. ENERGY-LEVEL SHIFTS NEAR A 
HALF-SPACE 

A. Ground state 



Substituting the photon propagator ([97]) into Eq. (|108|) 
we find that the energy shift of the atomic ground state 
\g) is given by 



AE, 



dq\\ q\\ 



duj 



9 yJq1+UJ^ 

(111) 



where we have used the notation 
and l/i^ip 



mil iMmil "'"iMmil 

MmiP- The reflection coefficients are as 
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defined in Eq. 
read 



In terms of the new variables they 



;tTE 



-TM 



(112) 



Note that we have replaced ^(w) by e{uj) in Eq. (|112p 
compared to Eq. ([88|) . because for the relevant frequen- 
cies both functions coincide (see Section IIVBI and Ap- 
pendix |B]). If we now introduce polar coordinates accord- 
ing to, Lo = uJmg X COS (/), = ujmg X sin (j) and then write 
y = cos 0, we obtain the perhaps most useful expression 
for the ground-state shift, especially for numerical anal- 
ysis and for investigating the effects of retardation, 



1 



dxx 



dy 



L0„ 



2(1 



1 -I- x^y^ 



(113) 



The result in Eq. ()113p formally takes the same form as 
the results obtained in calculations involving only non- 
dispersive dielectrics (see e.g. [l^), the only difference 
being the reflection coefficients that now, through the 
dielectric constant, depend on the product xy of the in- 
tegration variables which is the photon frequency in units 

of 



~TE ^ 1 - \J V'^W'^mgXy) - 1] + 1 

1 + ^ y'^{t(iu:mgxy) - 1] + 1 ' 



~TM ^ e(»^mga;y) - ^J y'^\t{iuj^gxy) - 1] + 1 ^^^^^ 
e{iuj„igxy) + y'^[e{iujrn.gxy) - 1] + 1 



Eq. ()113p is suitable for numerical analysis but does not 
give immediate insight into the dependence of the energy 
shift as a function of the distance from the surface. It 
is therefore instructive to consider some of its limiting 
cases. 

As has been spelled out e.g. in Ref. [l^, the dimension- 
less parameter that plays a decisive role in the charac- 
teristics of the Casimir-Polder interaction is given by the 
combination 2cj,„gZ which is the ratio of two time-scales: 

(i) the typical time 2Z/c needed by a virtual photon to 
make a round trip between the atom and the surface, and 

(ii) the typical time-scale a;"^ at which the atomic sys- 
tem evolves. While Eq. pi3|) includes a sum over atomic 
states |m), in reality contributions to the shift will be 
dominated by the state which is connected to the ground 
state by the strongest dipole transition. We shall call the 
frequency Wmg that pertains to this strongest transition 
the "typical transition frequency" and it is this number 
that enters the retardation criterion parameter. Roughly 
speaking, if 2LUmg2 <C 1 we are in the so-called nonrc- 
tardcd regime when the time needed by the photon to 



travel between the dielectric and the atom is negligibly 
small compared to the typical atomic time-scale. Then 
the interaction can safely be approximated as instanta- 
neous and our result should reduce to that calculated 
by Barton [IJl, who considers only the Coulomb interac- 
tion of an atom with surface polaritons. In the opposite 
case, 2uJmgZ ^ 1, the interaction becomes retarded, i.e. 
by the time the photon has completed a round trip, the 
atomic state has changed significantly. In that case, for 
reasons that are not obvious but will become apparent 
later, the interaction depends only on static polarizabil- 
ities, i.e. the polarizabilities evaluated at zero frequency. 
The diagonal polarizability of the spherically symmetric 
atom is then 



2E 



Wo 



and the susceptibility of the dielectric becomes 



(115) 



6(0) 



= 1 



as explained in Appendix IbI 



= 1 + -^, (116) 

Wrp 



Nonretarded limit 



In order to take the non-retarded limit of the energy 
shift it is best to start from Eq. (jllip . After changing 
variables from w to s with lo = (2w,„gZ(7|[) s we take the 
limit 2uJmgZ — > 0, which we may do because the line s = 
oo does not contribute to the integral, and approximate 



gjj [l + {2iU,ngZ)^S^] 

[l + e{2LUrngZfs^] « g|| . 



This significantly simplifies Eq. (|llip . The integral 
becomes elementary, and the final result reads 

e{iuj) — 1 
e{iuj) + 1 

(117) 



1 



327r2eoZ' 
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duj 



^mg 
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We observe the expected Z^^ behaviour of the energy 
level shift in the non-retarded or "van dcr Waals" regime. 
In order to see that the result in Eq. (|117p is equivalent to 
the slightly more awkward principal- value integral given 
in Eq. (7.14) of Ref. [H or Eq. (13) of Ref. [HI, one 
needs to re-write 



1 



1 



and the re-rotate the contour from lo to iu in the first and 
to —iuj in the second summand. Eq. (|117p also confirms 
the results derived on the basis of the phenomenological 
noise-current approach to quantum electrodynamics with 
dielectric media, see e.g. 
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2. Retarded limit 

In order to work out an approximation to the en- 
ergy shift when retardation is dominant, it is convenient 
to start with Eq. pi3p where the decisive parameter 
2uJmg2 is present in the exponential which in the Umit 
2ujrngZ — > oo strongly damps the integrand. Then the 
main contributions to the integral come from the neigh- 
bourhood of X = 0^ and one can obtain an asymptotic 
expansion of the integral by expanding the integrand in 
a Taylor series around this point. A straightforward cal- 
culation gives 



647r^eo 



E E 



u;2 Z5 



(118) 

where in the parentheses we have neglected terms of or- 
der Lo^gZ~^ and higher. The coefficients C45, to be 
given below, depend only on the static dielectric con- 
stant of the material. The fact that to leading order 
the Casimir-Polder force depends only on the static po- 
larizability of the atom, Eq. (|115p . is well known jSlj . 
Therefore the leading-order Z^^ term in Eq. (|118p is 
identical to the retarded limit of the energy shift of a 
ground-state atom interacting with a non-absorptive di- 
electric half-space described by a static refractive index 
n(0) = n = ^Jl + / , which has been derived previ- 
ously [1^. We just quote the results for the coefficients 
cP from Ref. [11: 
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n2 - 1 V 3 
2n4 



-n + n 
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Vr? + T+ 1 



{rfi - l)\/n2 + 1 In [\/n^ + l + n] j 
2n** - 2?i2 - 1 / /— N 



4 2 , 4 



In 



+ 1 + 1 



(n2 - 1)\A?TT \n [s/n'^ + 1 + ' 
2n2(2n4 - 2^2 + 1) 



In ( ri^ — 1 



(„2 _ 1)3/2 



In other words, to leading-order, in the retarded limit, ab- 
sorption makes no difference and only static polarizabili- 
ties, of both the dielectric and the atom, matter. This is 
because the photon wavelengths that matter the most in 
the atom- wall interaction are of the order of the distance 
between the atom and the surface of the dielectric and 
longer. Thus for an atom in the so-called far-zone only 
long wavelengths of the electromagnetic radiation come 
into play, which means low frequencies. 

Now we turn our attention to the next term in the 
asymptotic expansion which is proportional to . This 



is the first term that contains information about correc- 
tions to the energy shift due to absorption in the retarded 
regime. Apart from the factor 47/ w^, the dimensionless 
coefficients Cg depend again only on the static refractive 
index n = -^1 + Wp /w^ and are given by 



3(n- l)(7i+ l)2(n2 + 1) 



x<! 6n^ - 371^ - 1171"* 



-67T,^(7T,''' + 77'* - 77^ 

4 



477'' + 277^- 577 + 7 



277 - 2 ) In 



3(77- 1)(77+ 1)2(772 + 1) 



X <^ - 677* + 377^ + 1077^ - 577^ + 377"^ - 77^ - 677^ 



377 — 3 ) In 



+377-* (277^ + 277^ 



(119) 



We provide plots of these functions in Fig. [2] from where 
a quick estimate of the value of these coefficients can 

be obtained. Since both Cg and arc positive we see 
that absorption reduces the magnitude of the ground- 
state energy shift by an amount that is proportional to 
the damping constant 7, cf. Fig.O We also note that the 
correction goes with the inverse square of the absorption 
frequency cjt in the dielectric so that only absorption 
lines that lie at sufficiently low frequencies make a sig- 
nificant difference. This happens because the main con- 
tribution to the shift of the ground state in the retarded 
limit comes from long wavelengths or equivalently small 
values of x (which is a scaled frequency). Therefore, the 
integral is not sensitive to any absorption peaks which 
lie at higher frequencies as there the integrand is highly 
damped anyway, cf. Eq. pi3p . 



B. Excited states. 

The shift of an excited energy level gets contributions 
from both parts of A£'"\ Eqs. (fTOg]) and ^09]) . The non- 
residue contributions, Eq. (|108p . assume exactly the same 
form as the results of the previous section. Therefore we 
will not analyze them again but shall instead have a closer 
look at the additional contributions due to Eq. (|109p . 

Plugging in the photon propagator, Eq. ([97|) . we find 
that the energy shift of the excited state \i) is given by 
the real part of the following expression 
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(120) 
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FIG. 2. Plot of the coefficients (n) that enter Eq. (fTT8|) for 
different values of the static refractive index n = n(0). 



variable in Eq. (|120p to = xj^mi - qll\^mi\ and get 



87reo 



E 



X {(r™ - fc^-™) ImLP + 2 (1 - fc.^) ,(121) 

where the contour of integration runs from /c^ = 1 along 
the real axis to fc^ = and then up along the imaginary 
axis to kz = ioo. The reflection coefficients expressed as 
frmctions of k, are 



k, + ^[e(|w„„|) -l] + kf 



^TM^^ -) ^ 4\^,m\)kz ~ y^[e{\0J,m\) - 1] + fcg ^^^^^ 

e{\uj,m\)k, + V[e(|w,„„|) - 1] + fc2 ■ 

We now go on to analyse A£,f in the nonretarded and 
retarded limits. 



Nonretarded limit 
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FIG. 3. Plot of the exact ground-state energy shift (contri- 
butions due to the perpendicular component of the atomic 
dipole moment) ASg, Eq. multiplied by 2* ^mg as a 

function of Zujmg for various values of the damping param- 
eter 7. Solid line represents the energy shift caused by the 
non-absorptive and non-dispersive dielectric half-space with 
static refractive index n(0) — \/2. 



Here r^^ are the reflection coefficients of Eq. ([55]) evalu- 
ated at the atomic transition frequencies w = |w,ni|- Also, 
the restriction of the sum over atomic states to those ly- 
ing below the state \i) should be noted. For the purposes 
of asymptotic analysis of /S.8* we change the integration 



In the nonretarded limit of Eq. (I12ip we have 
2\ujmg\Z ^1. It is expedient to split the integration 
in Eq. (|12ip in the following way 



dk. 



d{ikz) 



dk,. 



(123) 



and note that in the limit 2\LiJrai\Z — > the second in- 
tegral on the RHS contributes to the asymptotic series 
only terms that are proportional to non-negative powers 
of Z and can therefore be discarded. The remaining part 
is given by 
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X ^(Cf + fc.'?™)l/^Lf + 2(l + fc.^)r™|/.;^J^] (124) 

where are the reflection coefficients of Eq. (|122p eval- 
uated at imaginary argument ~ f^^^^{ikz)- Scaling 
the integration variable according to x 
approximating 



-1]- 



(2|w„„|Z)2 2\uj„ 



(125) 



we derive that, in the nonretarded limit, Eq. ()120p be- 
comes 



as: 
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(126) 
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To leading order the residue contributions to the energy 
shift of the excited state cf. Eq. (|109p . are given by 
the real part of the above expression, 



AE7 



1 



E 



327reoZ3 |e(|wr„i|) 
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+ iF 



.(127) 



Thus in the nonretarded regime the residue contributions 
behave as and therefore are of the same order as the 
non- residue contributions, cf. Eq. (|117p . The result in 
Eq. ()127|) is in fact equivalent to the real part of Eq. (7.10) 
derived in Ref. 1341. 



2. Retarded limit 



the energy shift 



47reo \n{uJrni) + IP 
(|n(|a;,„i|)|^ - 1) cos(2|a;,„i|Z) 



-2Im[7i(|w,„i|)] sm{2\uJmi\Z) 
(|n(|w„„;|)p - l)sin(2|a;,„i|Z) 



-|-2Im[n(|a;mi|)] cos(2|wmi|Z) 



(2|c^™|Z)2 



(130) 



Now we turn our attention to the asymptotic behaviour 
of Eq. p2ip in the retarded limit, i.e. when 2|w,„i|Z ^ 1. 
It is again useful to split the integration in the same way 
as in Eq. (jl23|) . only that now both integrals play an im- 
portant role. The first contribution, the integral along 
kz S [0, ioo], given in Eq. (|124p . can be tackled by use of 
Watson's lemma [s^. Noting that for 2\uJmi\2 ^ 1 the 
integrand is strongly damped, we separate off the expo- 
nential and expand the remaining part into Taylor series 
about = 0. The resulting integrals are elementary and 
we obtain for the leading term 
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We see that in the retarded regime the two contribu- 
tions to the shift of an excited state behave quite differ- 
ently. The non-residue contribution in Eq. (|108p behaves 
as (see the analysis of the ground state shift in Sec- 
tion [VT^, and the residue contribution in Eq. (|130p de- 
pends on distance as . While it would be tempting to 
jump to the conclusion that Eq. (|130p will always dom- 
inate, this might in fact not always be the case as the 
relative size of the two contributions also depends on the 
values of the dipole matrix elements involved, which can 
vary significantly. Furthermore, Eq. ()130p is oscillatory, 
so that at least in principle there are sets of parameters 
for which it vanishes. Finally, we remark that it is easy to 
verify that in the limit of non- absorptive dielectric media 
our results reduce to those derived in Ref. [l9|. 



Next we deal with the integral on the interval k^ G [0, 1] 
which, unlike in the nonretarded case, cannot be dis- 
carded. However, its asymptotic expansion in inverse 
powers of Z is easily obtained by repeated integration by 
parts. Interestingly, the asymptotic series contain non- 
oscillatory terms that exactly cancel out the contribu- 
tions given in Eq. (|128p . Altogether we find that the 
leading-order and next-to-leading-order terms are 
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with the refractive index n(|wmi|) = \J e(|a;mi|). It is 
interesting to observe that to leading-order in Z only 
contributions due to the parallel component of the atomic 
dipole moment are contributing; contributions due to the 
perpendicular component of the atomic dipole moment 
appear only in next-to-leading order. Again we need to 
take the real part of Eq. (|129p to get the explicit form of 



VII. SPONTANEOUS DECAY RATES NEAR A 
HALF-SPACE 



The spontaneous decay rates are given by the imagi- 
nary part of the complex self-energy, Eq. (|110p . As the 
non-residue contributions to the self-energy in Eq. (|108p 
are real, these contribute towards the energy-level shifts 
only, and the decay rates are contained solely in the 
residue contributions to the self-energy, Eq. (|120p . which 
are complex. In the non-retarded limit the decay rates 
are given by the imaginary part of Eq. (jl26p . 
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FIG. 4. Normalized lifetime ry, Eq. (|133p . of the atomic state 
|i) plotted as a function of \ijJmi\/uJT- The sequence of graphs 
corresponds to various distances of the atom from the mirror 
Zujt, as indicated. The three different line styles (colours) 
indicate distinct choices of the damping constant 7 in the 
dielectric: 'j/ujt — 0.05 (black, solid), j/uit = 0.5 (blue, 
dashed), "//ujt ~ 5 (red, dot-dashed). For sufficiently high 
frequencies |cjmi| the dielectric becomes transparent. If the 
atom is close to the surface and absorption is small, the inter- 
action is resonant at |a;mi|/t»JT ~ 1, i-e. when the frequency 
of the atomic transition coincides with the absorption line of 
the dielectric. 



and in the retarded limit by the imaginary part of 
Eq. (fT29l) : 
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The result in Eq. p3ip is found to be in an agreement 
with that derived in Ref. [13, their Eq. (128). A con- 
sistency check on Eg. (|132p is that it reduces to the re- 
sults given in Ref. |19| if we assume n(a;) to be real and 
frequency- independent . 

As a numerical example we plot the normalized lifetime 
of the atomic state |z) decaying into a lower state |m) . For 
simplicity we assume a two- level system and l/i^J = 
so that the atom is polarized horizontally with respect to 
the surface. Then the normalized lifetime that we plot 
in Fig. m is given by 



AT, 
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where the quantity A£* comes from Eq. (|120p and AF? 
is the well-known decay rate in free-space 
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(134) 



VIII. SUMMARY AND CONCLUSIONS 

We have shown that starting from a gauge- independent 
microscopic model as represented by the Hamiltonian 
(fT9)) it is possible to develop a formalism which allows 
to calculate QED corrections in the presence of absorp- 
tive and dispersive boundaries. We have used a diagram- 
matic technique to integrate out the damped polaritons 
in order to arrive at a Dyson equation for the electromag- 
netic displacement field propagator. We have solved this 
integral equation exactly, using traceable methods. The 
knowledge of the exact propagator has enabled us to cal- 
culate analytically the one-loop self-energy diagram for 
an electron bound in an atom near a dielectric half-space 
and hence to determine its energy-level shifts and the 
change in transition rates, which derive from the real and 
imaginary part of the electron's self-energy, respectively. 
This is serving as a proof of principle that the theoretical 
framework developed here works correctly and efficiently, 
as most of these results have, in one form or another, been 
derived previously by other methods, though often with 
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considerably more effort or much less rigour, especially 
as regards basic principles. 

We have looked at the role of the material's absorp- 
tion in some detail and confirmed the previously known 
result that absorption has the most profound impact on 
the atomic system in the non-retarded regime, i.e. when 
the distance between the atom and the mirror is much 
smaller than the wavelength of the dominant atomic 
dipole transition. If the distance between the atom and 
the surface far exceeds the wavelength of this dominant 
transition, then, to leading order, dispersion and absorp- 
tion do not affect the ground-state shift for which only 
static polarizabilities matter. The next-to-leading order 
corrections are proportional to damping constant of the 
Lorentz-type dielectric function, and it turns out that 
only the material's absorption lines that lie in the low- 
frequency region have a significant impact on the ground- 
state energy- level shift. We have also rcderivcd the dis- 
tance dependence of the excited energy-level shifts and 
spontaneous decay rates. We have confirmed the fact 
that in the non-retarded regime the absorption is of fun- 
damental importance to both the change in decay rates 
and the energy-level shifts. For example, for an atom 
near a nondispersive dielectric the spontaneous decay 
rate in the near-zone comes out as a distance- independent 
constant [l9| , whereas in reality, when the absorption is 
taken into account, a distance dependence cx is ob- 
tained. In the far-zone or retarded limit, the presence 
of absorption does not affect the characteristic be- 
haviour of the excited energy-level shift and spontaneous 
decay rates, even though the coefhcients differ from the 
nondispersive case. 
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Appendix A: Dressed photon propagator as a 
boundary-value problem 

We aim to show that the integral equation satisfied by 
the dressed photon propagator derived in Section [IVBI 

D,,{ry;Lo)^D[l\v-v'-Lo) 

+ f d^r,gir,)D'^f (r - n; c.)i?,fc(ri, r'; uj). (Al) 

can also be solved by considering it as a boundary-value 
problem. Recah that D^l\r — r') is the photon propa- 
gator in free space, Eq. ([SU]) . and ^(r) is a dimensionless 
coupling constant that is equal to unity in the region oc- 
cupied by the dielectric and vanishes otherwise. To de- 
scribe a dielectric half-space occupying the z < region 
of space, as illustrated in Fig. [H we take ^(ri) = 9{—zi) 



where 6 is the Heaviside step function. Knowing that the 
free-space propagator satisfies the differential Eq. ([39]) we 
apply the same differential operator to Eq. (jAl[) . and af- 
ter a short calculation we obtain the differential equation 
satisfied by the photon propagator in the half-space ge- 
ometry. 
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The RHS can be re-writtcn as 



Djk{v,v'-uj) (A2) 
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where ^j^(r — r') is the transverse delta- function. Now 
it is more apparent that the RHS of Eq. (jA2p is a dis- 
tribution which, unlike the transverse delta function, is 
sharply localized around the point r = r', because the 
non-local part of the transverse delta function is removed 
by the application of the Laplacian as is obvious from the 
relation 



(A3) 



The locality of V^iSj^. (r — r') is very helpful towards the 



solution of the differential Eq. (|A2p . which is essentially 
a scattering problem. Its RHS contains a distribution 
representing a point-like source and our task is to work 
out reflection and transmission at the boundary of the 
dielectric. In order to proceed any further, we need to 
specify physical situation, i.e. decide on which side of the 
boundary the source is placed. Since our ultimate aim is 
to work out the energy-shift in an atom located outside 
the dielectric, we choose to consider the case z' > 0. 
Then we write Eq. (jA2[) in a piecewise manner; on the 
vacuum side we have 

-L.2^,fc(r,r';c^) = 6oV24(r-r'), z > 0, (A4) 
and on the dielectric side we have 

-a^)u:^D,k{v,v';uj) = Q, z<Q, (A5) 

with the behaviour of the propagator Djk{v,r'\oj) across 
the interface z = is still to be determined. The local 
character of the RHS of Eq. (|A2p simplifies its solution 
in that it makes the RHS of Eq. (jA5[) go to zero. 

In order to solve Eqs. (jA4[) and (|A5[) we start with the 
following ansatz 



Afc(r, r';w) = 



(r,r';c.) 



d'§{v~v'-uj) + DI'{v,v';uj) z>0. 



z < 0, 



{r)i 



(A6) 
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On the vacuum side we write the solution as a sum that 
consists of a particular solution ofl} (r — r'; w), which we 
already know from Section IlII Bl Eq. (|42p . and a solution 
D^^j) (r, r'; w) of the corresponding homogeneous equation 
(i.e. Eq. (|A4p with the RHS set to zero) which repre- 
sents the correction due to reflection at the boundary. 
The solution on the dielectric side Z?^^' (r, r'; w) repre- 
sents the transmitted part and satisfies the homogeneous 
Eq. (|A5p . The homogeneous solutions r'-]J^(r, r'; w) and 

o'f^ [t^t' ]uj) are chosen in such a way that the general 
solution in Eq. ()A6|) satisfies appropriate electromagnetic 
boundary conditions across the interface z = 0. To see 
what these boundary conditions should be recall the for- 
mal definition of the dressed propagator 



Eq. ((651) . nr^-y be written as 

-eo(q'-V2)5,^(q||,z-z'), z>Q (AlO) 
(vl-qjj+eH^') Aj(z,z') = 0, z<Q (All) 

where 8^^ (q|| , z — z') is the Fourier transform of 5fj (r — r' ) 
with respect to ru — ry 



5,:^(q||,z-z')- jS{v\\-v\)^ 



-•iq|r(r||-r' 



«)<5,^(r-r'). 



Therefore the homogeneous solutions D^^^ {v,v' ]uj) and 
D,j{v,v'-Lu) = -Un\T [D,{v,t)D,{v\t')] \^). (A7) ^fk^^^^''^^) Eq. dMl) must necessarily take the form 



fi 

The displacement operator Di{r,t) satisfies Maxwell's 
equations which follow from the Heisenberg equations of 
motion for the field operators. Therefore, the photon 
propagator, by virtue of its definition (jA7p . when taken 
as a function of argument r and index z, is required to 
satisfy Maxwell's boundary conditions across the inter- 
face: 



En continuous 



D, continuous 



Bii continuous 



D 



z=0- 



z=0+ 



2=0- 



z=0+ 



(A8) 

with 11= {x,y}. 

The apparent complication arising from the appear- 
ance of a non-standard distribution in the boundary- 
value problem (jA2[) is just an illusion. In fact, it is easier 
to find the solution of Eq. (jA2p than it is to solve the 
differential equation satisfied by the Green's function of 
the standard wave equation (see e.g. [l^l). Eqs. (|A4p and 
(|A5p . together with the boundary conditions (jASp . form 
a boundary-value problem which is equivalent to the in- 
tegral equation ()Aip with the choices g{ri) = 9{—zi) (di- 
electric occupying the left half-space) and z' > (source 
located in vacuum). 

In the following we shall use Eq. ([80]) in the process of 
matching the boundary conditions. This is safe because 
we consider z = 0='= and the source located at z' is always 
well away from the boundary, so that z z' is assured. 

To proceed further we note that taking the divergence 
of the integral equation in Eq. (|Aip and using the fact 

that the free-space propagator is transverse, ^iD^^ (r — 
r'; LLj) = 0, one infers that the dressed photon propagator 
is transverse everywhere as well 



V,A:fe(r,r';Lj) = 0. 



(A9) 



With this Eqs. (jA4p and (jASp simplify further and par- 
tially Fourier transformed into the (qy , z) space, cf. 



D'^iz.z') 



with kz 



^[i?.,(z')e^'=^^ + S.,(z')e-'^^^], z>0, 

(A12) 



[T„-(z')e-*'=-<^^ + U,j{z')e'''-''] , z < 0, 

(A13) 



qf. + IT] and k^d 



and the square roots taken such that Im(fcz) > and 
\m(kzd) > 0. With this choice of sign for the square 
roots, the terms in Eqs. (|A12p and (|A13p that contain 
exponentials e^*'"'=^ and e^^"^^- are unphysical as they rep- 
resent waves that diverge at infinity. Thus we must set 
Sij = = Uij. The remaining two matrices Rij and Tij 
are determined by the requirement that Eq. (|A6p satis- 
fies the boundary conditions in Eq. ()A8p . We note that, 
in addition, the transversality of the dressed propagator, 
Eq. (jA9p . imposes rather stringent constraints on both 
Rij and T^ . For example, the matrix Rij needs to be of 
the form 

Rtj =Vi{ci\\)rj{ci\\,z'), (A14) 
where the vector v is such that q • v = 0, leading to 



(A15) 



with q = (q|| , fcz). One might pick Vx ~ —Qy and Vy = 
so that 



{-qy,qx,0) . 



(A16) 



However, this choice is too restrictive on its own, as there 
is no a priori reason for Dzj to vanish. Therefore, an ad- 
ditional basis vector is needed in order to span the am- 
plitude Rij in full generality. An obvious and convenient 
choice is to choose a vector that is orthogonal to both q 
and V, 



w = V X q = (fefcz,gyA:^,-q|). 



(A17) 
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Then we can represent Rij as the hnear combination 

R,j = [av + /3w],,r,(z') 



(fc.)r™(z') + e™(fc.)r™(z')- 



where we have recognized, apart from normahzation fac- 
tors, the transverse electric and transverse magnetic po- 
larization vectors of Eq. (|78|) . Similarly we have 



T., = en~K,)tJ^iz') + (A18) 

We have chosen write out the fcz-dependence of the polar- 
ization vectors, even though and k^^ are expressible in 
terms of the frequency uj and the parallel wave- vector q|| , 
because this explicitly indicates the wave- vector to which 
a given polarization vector is orthogonal to. The de- 
composition into transverse electric and transverse mag- 
netic components significantly simplifies the matching of 
boundary conditions. The dressed photon propagator 
can now be written in the form 



k, 



9(z) 
(A19) 



the last term of which is the free-space photon propagator 
from Eq. ([80|) for z — z' < 0, as appropriate for the 
matching of boundary conditions at z = when z' > 0. 
Imposing the boundary conditions of Eq. (jASp wc find 
that 



Akzz' 



tx e^{-k,) 



with rx and tx being the standard Fresnel's reflection and 
transmission coefficients listed in Eq. ([55)1 . Plugging the 
above amplitudes into Eq. (|A19p we then readily obtain 
the photon propagator given in Eq. (j96p . 

Wc can readily apply the same methods to obtain the 
propagator in the case when the source is placed in the 
dielectric, i.e. for z' < 0. The calculation goes along 
exactly the same lines as for z' > and one can show 
that the photon propagator in the case of the source being 
placed in the dielectric is given by 



xlei-z) [^(c^)e,^(q||,-fc,,)e^^(q||,fc,,)ri] e-*'=-(^+^') + ^(z) [e,^(q|| , ft,)e^^(q|| , e^^.-^^..-' j, .(A20) 



Here the refiection and transmission coefficients are those 
appropriate for left-incident modes; they are given by 



tion relations (fn)) - p^ : 



„TE 



+TE 



„TM 



kz + k 
2k 



i{Lo)kz + kzd ' 



d ^TM ^ 2^g(a;)fc^d 
kz+k^d' ^ £.{uj)k, + k^d' 



It is easily verified that £'y(r, r';w) is indeed transverse 
everywhere. 



Appendix B: Simple model for t{Lj) 

In order to determine the dielectric permittivity of our 
model we use the equations of motion for the fields that 
follow from the Hamiltonian ([7|)- ([TT|) and the commuta- 



d_ 
dt 
d_ 

dt 
d_ 

dt 
d_ 
dt 



d 



D(r,i) 
B(r,t) 
X(r,t) 
P(r,t) 



Mo 



-V X B(r,t) 



1 



-V X E(r,t) 

-Ala;|X(r,t) +g(r)E(r,t) 



/ p,V^Y,{Y,t) 

Jo 



1 



-Y,(r,t) = — Z,(r,t) 



d 



(Bl) 
(B2) 
(B3) 

(B4) 
(B5) 



^Z,(r, t) = -p,iy^Y,{r, t) + p,i/2X(r, t) (B6) 
ot 

First we deal with the subsystem consisting of the po- 
larization field and the reservoir. It is well known |24l | 
that when a quantized harmonic oscillator is coupled to 
a bath its equation of motion takes the form of a quan- 
tum Langevin equation. Thus we expect the equation of 
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motion for the polarization field, which is nothing but a 
set of independent oscillators, to take the form 

+ MulX{r, t) = F..an(r, t) + Fext(r, t) 

(B7) 

where iJ,{t—t') is the so-called memory function related to 
dissipation and Fy„„(r, t) represents some random force 
operator (see e.g. [23] for details). Both ii{t — t') and 
Fran(r,i) arisc as a consequence of the coupling to the 
bath and are to be determined in terms of the parameters 
of our model. The term Fgxt {t) represents any external 
forces (i.e. those in addition to the harmonic restoring 
force) that may be applied to the polarization field. To 
show that Eqs. (|B3p - (|B6p indeed combine to yield an 
equation of the form of Eq. (jB7|) . we eliminate P(r,t) 
and Z^(r,t) and rewrite the equations for X(r,t) and 
Yy(r, t) as 

(^^+z.2)Y,(r,t) = ^2X(r,t). (B9) 
The most general solution of Eq. (jB9p may be written as 

/oo 
dt'G,it^t')X{r,t) (BIO) 
-oo 

where Y^(r,<) is the solution of the homogeneous equa- 
tion, i.e. Eq. (|B9|1 with its RHS set to zero. 



(r,t) = Y.(r,0)cos(.i) + Ml^El^lM. (bh) 

We assume that bath operators Y,y(r, 0) and Zi/(r, 0) sat- 
isfy the canonical commutation relations at the initial 
time t = 0, ci. Eq. (|15p . which we take as a moment in 
the distant past when the interaction has been switched 
on. The second term in Eq. (jBlOp is a particular solution 
expressed in terms of the Green's function 



did 



z/^ — (oj + iey 



j(t-t') 



(B12) 



The ie prescription for handling the pole ensures that we 
have a retarded Green's function with Gi,{t — t') = for 
t-t' <0. Eq. pT2|l is easily obtained from Eq. (|B9| by 
using Fourier transforms according to 



Y,(r,c^) 



dte*"*Yj,(r,t). 



(B13) 



Note that the choice of the retarded solution breaks time 
reversal invariance, as has been noted in Ref. [l^l . The in- 
tegral in Eq. (jB12p is easily worked out using the residue 



theorem, and Eq. (jBlOp may be rewritten as 
Y,(r,t)=Yf(r,t)+X(r,t) 

- [ dt' cos[iy{t-t')]^JL{r,t) (B14) 

where we have integrated by parts. Plugging the above 
expression into Eq. (|B8P we obtain 

Al^X(M) 

+ J d<' 1^ di^p,,y^ cos[,y{t - t')] I ^X(r, t') 

poo 

+7Wtj|X(r, t) = / di^p^,y^Yl^{r,t). {B15) 
Jo 

This is the quantum Langevin equation that follows from 
our model. Comparing with Eq. (|B7p lets us to identify 



Fran(r,t) 
M(t - t') 



dvp,v'Y^{v,t) 



(B16) 



v'^ cos[v {t-t')]. (B17) 



Now we are in the position to choose the bath oscillator 
masses pi,. Having in mind a simple single- resonance 
model of the dielectric permittivity, we choose pi, in such 
a way that the friction term in Eq. (jBlSp is local in time 
i.e. it is non- vanishing only for t = t' . This is achieved 
by choosing 



(B18) 



which gives a frequency-independent coupling between 
the bath and polarization oscillators, cf. Eq. ([TO]) . Then 
Eq. (jBlsp becomes 

— X(r, t) + 27^X(r, t) + c.^ X(r, t) 

= ^ j^dup,u^Y^{v, t) + ^E(r, i). (B19) 

We have augmented this equation by the "external force" 
term that arises when the polarization field is coupled to 
the electromagnetic field, which according to Eqs. (|Bip - 
(|B2p satisfies the equation of motion 



92 

V X [V X E(r, t)] + Aioeo^E(r, t) = -/io5(r)^X(r, t) 

(B20) 

with D(r,t) ^ eoE(r,i) + g{r)X{r,t). Similarly to the 
reservoir field discussed before, the most general solution 
of Eq. (|B19p is given as a sum of the homogeneous solu- 
tion (i.e. the solution of Eq. (jB19P with the RHS set to 
zero and the assumption that the oscillators are under- 
damped) and the particular solution. The homogeneous 
solution is of the same form as Eq. (jBlip except for an 
additional damping factor proportional to e~'^*. Since we 
assume that the initial time is a moment in the distant 
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past we may discard the homogeneous solution which is 
exponentiaUy smah for 7t ^ 1. The particular solution 
is easily obtained in Fourier space and is given by 



1 



2ttM 



'^^. F,an(r,a;)+g(r)E(r,a.) ^_^^, 



(B21) 

where Fi.an(i",w) is the Fourier transform of Eq. (|B16[) 
and is given explicitly by 

rra„(r, |w|) = 47Y|^|(r,0) + *-^|a;|Z|^|(r, 0). (B22) 



Substitution of the solution (|B2ip into Eq. (|B20p yields 
V X [V X E(r,w)] 

, gHr) 1 



9{r) 



Mo- 



fan(l",I^) 



E(r,L^) 

(B23) 



We may now read off the dielectric function given by 



- = l+.g2(r). 



2*70; 



(B24) 



with Wp = l/Aleo- The quantity that appears on the 
RHS is proportional to the so called noise-current oper- 
ator which is introduced ad hoc in the phenomenological 
formulation of the quantum theory developed in Ref. Q . 
In fact we have 



JN(r,w) 



. g(r) F,a„(r,w) 

-VjJ- 



(B25) 



Since the operator ri.an(r,'^) depends only on the initial 
coordinates and momenta of the bath, cf. Eq. (|B22p . for 
which the commutation relations are known, it is rela- 
tively easy to verify that 

Ji(r,cj), jt(r',a;')j = 'iT:ntolm[e{v,uj)]uj'^5^^\v - v') 

xS{lu -Lu')5,k . (B26) 

This derivation justifies these phenomenologically intro- 
duced commutation rules on a microscopic level. We just 
note that this result differs from Eq. (|C2ll to be used in 
the following Section by a factor of (27r)^ due to a differ- 
ent definition of the Fourier transform. 



Appendix C: Photon propagator from 
phenomenological QED 

The phenomenological theory of quantum electrody- 
namics, as developed in Ref. gives the electric field 
operator as 

E,{r,t) = -ino jd'v' j dwe^'"*G,fe(r,r';cj)Jfc(r',L^)+H, 

(CI) 



where Jk{r,uj) is the so-called noise current operator sat- 
isfying the following commutation relation 



TT 

x6iLj - Lj')5,k , (C2) 



and Gik{v, r';uj) is the Green's function of the wave equa- 
tion satisfying 

^S,kS^'Hr~r'), (C3) 

with the additional requirement that it is retarded in 
time. Note, however, that there is no transversality con- 
dition imposed and the RHS of Eq. (|C3P is just a diagonal 
S function. For an overview of the noise-current approach 
and some applications see Ref. In the following we 
shall use two properties of the Green's tensor in particu- 
lar, its reciprocity 



G,k{r,r'■,uJ)=Gk^{r',r■,iu) , (C4) 
and the integral relation 

fdV Im[e(r,c.)]G*,(r,r";c.)G,fc(r',r";c.) 

= Im[G,,(r,r';c^)].(C5) 



To prove the latter one multiplies Eq. ()C3P from the left 
by G^j(r",r;aj) and integrates over r. Then taking the 
difference between the resulting relation and its complex 
conjugate integrated by parts yields Eq. (|C5|) . 

In order to calculate the Feynman propagator of the 
electric field operator, i.e. the quantity 



Df^{r,r',t,t') = --{Q\T[E,{r,t)E,{r\t')]\0), (C6) 



we substitute into the above definition the operator (|C1 
and use Eqs. (IC2|-(IC5|. We arrive at 



dujoj^ \9{t ~ t')e-"^^*-'"> 
Ti'eo Jo L 

+e{t' -t)e"^^*-''A Im[Gy(r,r';a;)]. 



Now we carry out the Fourier transform with respect to 
t — t' using the distributional identities 



(C7) 



where V denotes the Cauchy principal value, and obtain 



9 foo 3 

Df^ir,r';n)^—V / du;—^ 2lm[G,,(r, r'; l.)] 



C. 



/ duju^[S{n-oj) + S{n + oj)]lm[G^j{r,r';u)lC8) 
eo Jo 

The Green's tensor must satisfy retarded boundary con- 
ditions in time in order to preserve causality. This means 
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that it is analytical in the upper half of the complex lo 
plane. Analyticity in the upper-half of the w-plane leads 
to Kramers-Kronig relations [1^, so that the Green's ten- 
sor inherits the causality properties of the permittivity. 
In particular, its imaginary part is an odd function of fre- 
quency a;, whereas its real part is even in u. With that 
we can proceed to deal with the principal- value integral 
inEq. Since Im[Gy (r, r'; -w)] = -Im[Gy (r, r'; w)] 

and the remaining part of the integrand is also odd, we 
extend the lower integration limit to — cx) and compen- 
sate by multiplying by 1/2. On the other hand, the real 
part of the Green's tensor is even in w, so that we can 
replace 



for whose calculation we have used the fact that asymp- 
totically the Green's tensor behaves as [11] 



lim Lo''G^j{v,v-uj) = -%,5(3)(r - r'). 



(CM) 



The 5 function integral in Eq. (jCSp is easily seen to be 



-O^ Im[Gy(r,r';|l]|)] 
Co 



(C15) 



so that the final result for the relation between the photon 
propagator and the Green's function of the wave equation 
on the real f2-axis can be compactly written as 



1 



Im[G,,(r,r',L^)] ^ -G,, (r, r', l^) (C9) 



without changing the value of the integral. Thus the 
principal-value integral in Eq. (jC8[) becomes 



T 



du- 



rG,y(r,r';w). 



(CIO) 



To work out this integral we consider a contour of inte- 
gration 7 that runs from — oo to oo and above the poles 
at w = ±n and then closes up in the upper half of the 
w-plane along the large semicircle |a;| — >■ oo. Because the 
Green's tensor in analytic in the upper half-plane, the 
such calculated integral vanishes and we can express the 
principal-value integral as 



V 



(Cll) 



where 7^ denotes the clockwise contours that go around 
the poles at w = ±$7 respectively and F denotes the 
contribution from the large semicircle taken counter- 
clockwise. Using the residue theorem we derive that the 
contribution from 7^ is given by 



Z?g.(r,r';0) = -— G.,(r, r'; |17|) - -S,^S<-^\r - v'). 

A similar formula has been given in Ref. [401 . We would 
like to use this result for a comparison with the results of 
Section flVBI First we need to emphasize that what we 
have calculated here is the propagator for the electric field 
E, whereas Section IIVBI derives the propagator for the 
displacement field D. Therefore, the results can coincide 
only when r and r' are both located outside the dielectric, 
which is why we restrict ourselves to this case. Then the 
Green's tensor Gij(r,r';w) splits into a free-space part 

gI^^ and a correction G ■'^^ that describes the reflection of 
the electromagnetic field from the surface, and Eq. (|C16|) 
can be rewritten as 



Gf{v-v'-\n\) + 5,,5(^\v-v') 



-Gi'')(r,r';|0|). (C17) 



This makes clear that the Feynman propagator is an 
even function of f2, unlike the Green's function of the 
wave equation which has the same analytical structure 
as the dielectric function. It is not difficult to verify that 
for the particular geometry considered here, the dielec- 
tric half-space, Eq. (|C17|) indeed holds. The terms in 
square brackets combine to deliver the transverse free- 
space propagator as given in Eq. The reflected part 



G|;^(r,r';|17|), which can be found e.g. in Ref. 



sat- 



isfies the homogeneous wave equation. Therefore, it is 
automatically transverse 



-n'G'.^ivy-n), 

eo 



(C12) 



The large semicircle F contributes the delta function 



i5.,5(3)(r-r'), 
eo 



(C13) 



V.G«(r 



r» = G 



(C18) 



and for real lo it coincides with the reflected part of the 
photon propagator Df- (r, r'; w) given in Eq. ([96|) . though 
away from the real axis they are different due to the 
different boundary conditions in time. Df^ir, r'; t — t') is 
a Feynman propagator whereas G^ (r, r' ,t — t') gives the 
retarded solutions of the wave equation. 
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